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Abstract 

Continuous spectroscopic measurements of a strongly driven superconduct- 
ing qubit by means of a high-quality tank circuit (a linear detector) are under 
study. Output functions of the detector, namely, a spectrum of voltage fluc- 
tuations and an impedance, are expressed in terms of the qubit spectrum and 
magnetic susceptibility. The nonequilibrium spectrum of the current fluctu- 
ations in the qubit loop and the linear response function of the driven qubit 
coupled to a heat bath are calculated with Bloch-Redfield and rotating wave 
approximations. Backaction effects of the qubit on the tank and the tank on 
the qubit are analyzed quantitatively. We show that the voltage spectrum of 
the tank provides detailed information about a frequency and a decay rate 
J> , of Rabi oscillations in the qubit. It is found that both an efficiency of spec- 

troscopic measurement and measurement-induced decoherence of the qubit 
demonstrate a resonant behaviour as the Rabi frequency approaches the reso- 
nant frequency of the tank. We determine conditions when the spectroscopic 
■ observation of the Rabi oscillations in the flux qubit with the tank circuit can 

be considered as a weak continuous quantum measurement. 
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I. INTRODUCTION 



Quantum electrical engineering treats electrical circuits as macroscopic quantum systems. 
Some of these systems are approximately characterized by two states, and, therefore, can be 
considered as prototypes of a quantum bit, a main element of quantum computers. Recently, 
an existence of two quantum states has been proven experimentally in electrical circuits based 
on Josephson junctions [1-7]; in so doing coherent oscillations between macroscopically 
different states have been measured both in a free-evolution regime (quantum beatings) 
and in the presence of resonant driving field (Rabi oscillations). In the majority of the 
experiments a detector, for example, a dc-SQUID, was strongly coupled to the qubit during 
the measurement that resulted in a fast collapse of the circuit into one of its eigenstates. For 
mapping a whole evolution of the qubit the measurements were repeated thousands times 
with the same conditions to gather a statistical ensemble. An alternative procedure has been 
employed in an experiment [7] where continuous measurements of Rabi oscillations in a three- 
Josephson-j unction (3 J J) flux qubit have been performed. In this experiment the qubit is 
inductively coupled to a high-quality tank circuit that serves as a linear detector. Rabi 
oscillations reveal themselves in a spectrum of voltage fluctuations in the tank as the Rabi 
frequency of the qubit Qr passes through the resonant frequency of the tank ojt- It should 
be noted that the approach undertaken in Ref. [7] can be considered as an experimental 
realization of weak continuous measurements that have been studied theoretically in Refs. 
[8-13]. Previously an effect of which-path detector on electron dephasing in a double-path 
interferometer was demonstrated experimentally in Refs. [14,15]. Weak coupling between 
the qubit and the detector (tank) does not yield complete information about the state of 
the qubit before and during the measurement while not introducing severe decoherence into 
the quantum system. Because of the last reason the qubit can be monitored continuously 
with an extraction of useful spectroscopic information about quantum processes going on in 
the system. Among other things, a Rabi frequency (VLr/2t{ = 6.284MHz) and a life-time of 
Rabi oscillations (tr^ = 2.5/xs) in the strongly driven qubit have been measured in Ref. [7]. 
We interpret the qubit as a strongly driven system if the frequency of Rabi oscillations, Or, 
exceeds significantly the damping rate of the qubit, Y : Q R ^> T, while being much less than 
the frequency of quantum beatings. Here we have to mention some distinctions between 
the theoretical analysis of the weak quantum measurements given in Refs. [8,9] and the 
experimental implementation [7]. An assumption made in the theory that a characteristic 
time of the detector is much shorter than the period of measured oscillations is not valid for 
the tank (detector) with the resonant frequency ojt of order of the Rabi frequency Qr which 
is measured in the experiment. This discrepancy is cleared up with a replacement of a local- 
in-time response coefficient of the detector by a nonlocal response function of the tank. In 
the output spectrum we will have now a product of the qubit spectrum and a Lorentzian that 
describes a susceptibility of the tank (a Fourier transform of the tank response function). 
The next distinction is related to the fact that the analysis performed in Refs. [8,9] is 
aimed at the continuous measurements of quantum beatings which are determined by the 
equilibrium spectrum of the qubit. This spectrum has a peak at the frequency of quantum 
oscillations between degenerate states of the qubit. By contrast, the experimental setup 
deployed in Ref. [7] is dealing with a strongly nonequilibrium situation when the measured 
quantum system is driven by the external microwave field. In this case the output spectrum 
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of voltage fluctuations in the tank (detector) depends on the nonequilibrium current (flux) 
noise in the qubit loop and has a peak at the Rabi frequency which is much lower than the 
frequency of quantum beatings but much higher than the decay rate of the qubit. Emission 
and absorption spectra of the driven atoms (two-level systems) have been investigated by 
Mollow [16] in the Markov approximation. It was shown that in the presence of a strong 
driving field a single line-shape function of the atom bifurcates into two Lorentzians shifted 
from the usual resonant frequency of the atom by the Rabi frequency which is proportional 
to the amplitude of the driving field. The signals from the atom can be detected by means 
of high-frequency spectroscopic devices only. The two-level system measured in Ref. [7] 
produces also a low-frequency output that can be associated with biasing of the qubit from 
the degeneracy point [17]. We suppose as well that relaxation and decoherence rates can 
not be introduced into the qubit equations phenomenologically as it was done by Mollow 
for the atomic system. The results of the experiment [6] are indicative of a pronounced 
dependence of the qubit damping rates on the amplitude of the driving force. In particular, a 
decay time of Rabi oscillations, r Rabi ~ 150ns, measured in Ref. [6], is significantly different 
both from a dephasing time {t 9 = 20ns) and from a relaxation time (r re i ax = 900ns) of 
the undriven flux qubit. It is shown theoretically [18] that well-known formulas for the 
dephasing and relaxation times [19,20] are no longer valid for the strongly driven qubit. In 
this case relaxation and dephasing are mixed and determined by the spectrum of the heat 
bath fluctuations S(cu) taken at the Rabi frequency as well as at combinations of the energy 
splitting of the qubit, co c = VA 2 + e 2 , and the Rabi frequency cu c ± £l R . Here A and e 
are a tunneling rate of the qubit and a bias, respectively. Besides that, the driven qubit is 
no longer in thermodynamic equilibrium with the heat bath as evidenced by zero value of a 
steady-state population difference between the qubit energy levels [18]. Because of this, we 
can not resort to the Callen-Welton fluctuation-dissipation theorem [21] to find the spectrum 
of qubit fluctuations. 

This paper is devoted to a detailed consideration of nonequilibrium fluctuations and 
decoherence in a strongly driven qubit coupled to a linear detector (a high-quality tank). To 
accomplish these ends we apply a formalism of quantum stochastic equations proposed and 
developed in Refs. [22-24]. Our quantitative analysis is motivated by the recent experiments 
[7,6] and based on an assumption of a weak interaction between the qubit and the detector. 
This assumption can fail near the point of exact resonance between the qubit oscillations 
and electromagnetic oscillations in the tank, in particular, at the point = u>t- We define 
conditions whereby the measurements of Rabi oscillations in the flux qubit performed in 
Ref. [7] with a high-quality LC circuit (a tank) fall into the category of weak quantum 
measurements [8,9]. With this aim in mind we calculate a spectrum of voltage fluctuations 
in the tank as well as a contribution of the detector into decoherence rate of the qubit. The 
experimental set-up implemented in Ref. [7] (so called impedance measurement technique) 
can also monitor an effective impedance of the system " qubit +tank" by applying a small 
ac current, hias, to the tank with a subsequent measurement of an angle between the tank 
voltage and the ac current [25,26]. This angle is determined by the magnetic susceptibility 
of the qubit provided that the frequency of the ac current coincides with the resonant 
frequency of the tank. Here we calculate the magnetic susceptibility of the strongly driven 
qubit together with its decay rates taking into account detuning between the high-frequency 
source and the qubit. 
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The present paper is organized as follows. Dynamics and fluctuations in the linear 
detector (a tank circuit) coupled to the qubit are under study in Sec. II. We derive expressions 
for the averaged voltage in the tank, for the angle between the voltage and ac current, as 
well as for the spectrum of voltage fluctuations in terms of the magnetic susceptibility of 
the qubit and the spectrum of qubit fluctuations. In Sec. Ill we derive Heisenberg-Langevin 
equations for the driven qubit interacting with the tank and with its internal heat bath 
which is responsible for the qubit decoherence in the absence of the tank. These equations 
are subsequently simplified using a rotating wave approximation for the qubit that is weakly 
coupled to its environment, i.e. to the tank and to the internal bath. The decay rates of 
the qubit depending on the amplitude of the driving force and on its detuning from energy 
splitting of the qubit are derived in Sec. III. A dissipative evolution of a probability to find 
the qubit in the excited state as well as an evolution of an averaged current in the qubit loop 
are considered in Sec. IV. In the same section linear response functions of the driven two- 
level system are calculated together with parameters which are required for the impedance 
measurement technique (IMT) . The nonequilibrium spectrum of qubit fluctuations are found 
in Sec.V for the case of zero detuning. The output of the linear detector, namely, the 
spectrum of voltage fluctuations in the tank, and the contribution of the measuring device 
into qubit decoherence are presented also in Sec.V. In Appendices we outline our approach 
to the theory of open quantum system (Appendix A) and explain in more detail a derivation 
of collision terms (Appendix B) and spectra of fluctuation forces (Appendix C) . 



II. LC-CIRCUIT INDUCTIVELY COUPLED TO THE QUBIT 

We consider a 3JJ flux qubit [27,28] driven by a strong high-frequency field and induc- 
tively coupled to a tank circuit. This coupling is proportional to the coefficient of qubit-tank 
mutual inductance k^J ' L q L T as well as to the product of currents in the tank I T and in the 
qubit loop I q . Here k is a dimensionless coupling parameter, L q and L T are the inductances 
of the qubit loop and the tank, respectively, Ct is a tank capacitance. In the quantum case 
an operator of the qubit current is determined by a z matrix, I q = I q a z , that corresponds 
to two directions of the persistent current in the qubit loop. An operator of the current in 
the LC-circuit (tank) is defined in terms of creation-annihilation operators of photons in the 

tank having a resonance frequency uj t = 1/V L t Ct '■ It = /2L T (a + a + ), [a, a + ]_ = 1. 

For the Hamiltonian of the total "qubit-tank" system we have the expression 

A e 

H = —cr x + -a z - a z Fcosu t - a z (Q + / + \I T ) + H T + H qB (1) 



where A is a tunneling rate of the qubit, e is a bias, A = kI qS JL q L T is the coupling coefficient 
between the qubit and the tank. An operator Q describes an internal dissipative environ- 
ment of the qubit (without the tank), H q s is a Hamiltonian of this heat bath. The heat 
bath Q corresponds to all sources of external flux noise which are additional to the noise 
created by the tank [29,30]. In particular, the flux qubit is coupled to nuclear and impurity 
spins that can contribute to its decoherence and dephasing [31,32]. We also introduce here 
fit), a small external force that is required for a calculation of magnetic susceptibility and, 
in particular, an absorption spectrum of the qubit. In its turn, the tank driven by a bias 
current Iu as is characterized by the Hamiltonian H T , 
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H T = hu T (a + a + 1/2) - (a + a + )Q b - L T I T I Uas + H TB (2) 

Q b is a variable of another heat bath which directly interacts with the tank. That heat bath, 
having a free Hamiltonian Htb, is responsible for finite life time of the photons, t^ 1 , as well 
as for a finite quality factor of the tank, Qt = ujt/^t- 

Operators of the tank current, It, and tank voltage, Vt = i\JhuT/2CT(a + — a), obey the 

equations: It = Vt/Lt, 



From here on we put H — 1, and ks — 1. For a small coupling between the tank and its own 
bath and/or for the case of Gaussian fluctuations of free variables of this heat bath, Qf\ a 
total operator Qb(t) has the form 

Q b (t) = g[ 0) (t) + pL T /u T J dt^^t^lAh), (4) 

with 

Vb {t,t l ) = {i[Qf\t),Qt\t 1 )]_)d{t-t 1 ) (5) 
being a linear response function of the bath. Here 6{r) is the Heaviside step function, and 

(Qi 0) ) = o. 

To characterize this heat bath thoroughly we introduce also a correlation function of the 
unperturbed variables Q^: 

m 6 (mo = ((i/2)[gf ) (t),gf ) (ti)] + > (6) 

together with a corresponding spectral funcion Sf,(u) which represents a Fourier transform of 
M h {r). In the case of Ohmic dissipation in the tank with a resistance Rt the imaginary part 
of the susceptibility Xb(w), corresponding to the response function tpb(t — ti), is proportional 
to the frequency uo: x'li 00 ) = {1t/2^t)u, and to the line width of the tank 77- = 1/(RtCt) 
with v?&(t) = —(jT/^T)(d/dr)5(r). According to the fluctuation-dissipation theorem for 
the spectrum Sf,(u) we obtain: 

S b (u) = (7 T /2o; r )a;coth(a;/2T), (7) 

where T is an equilibrium temperature of the heat bath coupled to the tank. We suppose 
that this initial temperature is equal to the equilibrium temperature of the bath interacting 
with the qubit. 

In the presence of a time-dependent bias current, Ibias(t), the equation (3) for the operator 
of the tank voltage can be rewritten in the form of a simple stochastic equation 



l^ + lT 7t + u *) ^ T = V TT^ + Xuj2t&z + ~k hias - (8) 



T 
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A qubit operator a z in the right-hand side of Eq.(8) functionally depends on the tank voltage 
Vt- For a small qub it-tank interaction this dependence is approximately described by the 
formula: 



(9) 



Here we take into account the relation (5/ 5 It) = X(5/5f), as well as the equation: I T = 
Vt/ Lt- Fluctuations of the term & Zj0 , (<r z ,o) — 0, are determined by the internal bath of the 
qubit, Qo, only, so that the operator o z $ has no correlations with the heat bath Qb coupled 
directly to the tank. The functional derivative, (5a z (t)/5f(t')), has a magnetic susceptibility 
of the qubit, Xzz(w), as its Fourier transform: 



W) 



> 



duj 
2^ 



(10) 



These functions describe a behaviour of the qubit current, (I q ) = I q {cr z ), induced by 
var iations o f the time-dependent external flux, which can be created by the tank, / = 
k^jL q I 2 /L T & r f(t)- Taking into account a qubit back-action on the tank we obtain the fol- 
lowing equation for the tank voltage: 



d 2 



d 



dT 2+lT Jt +UJT , 



X 2 2/ Sa z (t) 



V T (ti) = 



-7~Qb > ^ + ^T&zfi + —hias- 
L/t 



1 



(11) 



The voltage Vr(t), or, more precisely, its average value and correlation functions, serves as 
a meter in the process of qubit measurements. Averaging of Eq.(ll) over the initial state 
of the qubit and over the equilibrium fluctuations of all dissipative environments allows us 
to find the average tank voltage, (Vr(t)) = Vycos(o;t + 0), induced by the time-dependent 
bias current, Ibi as {t) = I ac cosut. In the framework of the impedance measurement technique 
(IMT) [25,26] an imaginary part of the total impedance of the system " qubit +tank" , defined 
by a voltage-current phase shift 6, is studied as a function of qubit parameters, such as a 
bias, etc.. Matching the coefficients before e~ tuJt in the averaged Eq.(ll) gives the result: 



V T e~* e = 



-IUJ < UOn 



X 2 

i - y-xL' 

LiT 



It + 



X 2 uOj> a 



jr- (12) 



Here xL( u; ) anc l x'Li 00 ) are ^ ne Tea ^ anc ^ imaginary parts of the qubit magnetic susceptibility 
(10). At resonant conditions when the frequency of the bias current is exactly equal to 
the resonance frequency of the tank, uj = ut, the amplitude of the voltage oscillations is 
determined by the equation 



Vt = ^ {[k 2 L q I 2 q UTx' zz Mf + [It + k*L q lluTx" zz M?\ 



1/2 



(13) 



For the voltage-current phase shift we obtain the expression 



6 



tanfi __^ 2 u; TX ' zz M 



(14) 



Measurements of the angle between the average voltage in the tank and the bias current 
gives us an immediate information about the real part of the qubit magnetic susceptibility 
XzzO^t) taken at the resonance frequency of the tank ujt- A sensitivity of the impedance 
measurement technique is adversely affected by a qubit contribution to the tank damping 
rate which is proportional to the imaginary part of the qubit susceptibility x'L( u t)- 

It follows from the stochastic part of Eq.(ll) that a correlator of the voltage fluctuations 
(a cumulant function), 



M v (t,f) = ((l/2)[V T (t),V T (t')} 
satisfies the equation 

dti 



Z7T 



(15) 



d 2 d 9 \ A 2 9 ,8a z (t). 



dt 2 



dt 



I 



dto 



d 



d 



X 2 



^T 2 +lT^ + ^)5(t'-t 2 )-—u; 2 T ( 



dt 



W2) 



x 



d 2 f 2cu T 



The total spectrum of the qubit, S zz (u) = S zz ${uj) + S ZZ: t(uj), 

S Z M = J d(t-lf)e^*H±[a,(t),a s (lf)] + ^ 



(16) 



(17) 



incorporates a Fourier transform, S ZZt o(u), of the correlator (1/2 [cr Z fi(t),a Zt o(t')] + ), which is 
originated from qubit coupling to its internal heat bath, Q , together with a contribution 
S zz ,t(u), resulting from the qubit-tank interaction. This contribution has been built into 
the left-hand side of Eq.(16). 

It follows from Eq.(16) that the total spectrum of voltage fluctuations in the tank Sv(oo) 
(15) contains a contribution of the thermal noise Sb(u) (7) which is complemented by 
nonequilibrium noise generated by the qubit S ZZt o(u>) : 



S v {u) 



uj 2 — x 



2S b (iu) + k 2 LJ 2 u T S zzfl (uj) 



-,2 



[LO. 



UJ 2 ) 2 + UJ 2 ^t 



(18) 



Here ojt is a resonant frequency of the tank shifted in the presence of the qubit, 

u T = uj T ^l-k 2 L q I 2 qX ' zz M, (19) 

and 7t is a tank linewidth having regard to the qubit contribution to the tank damping, 

7t = 7T + k 2 L q I 2 uj TX ' zz {ujt) (20) 

A frequency shift of the tank and a correction to the tank damping rate depend on the 
total susceptibility of the qubit, Xzz(w), that should be calculated with consideration for all 
mechanisms of qubit dissipation. 
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III. QUANTUM LANGEVIN EQUATIONS. 



In this section we derive Heisenberg-Langevin equations with a subsequent goal of finding 
the nonequilibrium spectrum S zz {uj) of the qubit together with its magnetic susceptibility 
Xzziyi). To do that we consider a two-state system (a quantum bit) interacting with a 
heat bath Q in the presence of a harmonic driving force F(t) = F cos uj t. This heat bath 
incorporates a contribution of the internal qubit bath, Q , as well as a contribution of 
current (flux) fluctuations in the tank, XI T : Q = Q + XI T . An interaction with this bath, 
Hint = —Q&Z, has been integrated into the Hamiltonian (1). 

We suppose that the frequency of the external field uo can be different from the energy 
splitting of the qubit co c = VA 2 + e 2 with small detuning 5 = ujq—uj c , 5 <C uo c . In the rotating 
frame of reference the qubit is described by the operators: 

v A £ 
X = — a x H a z , 

f A £ \ ■ 
Y = <Xy cosuj t + — o z a x sin u Q t, 

/A e \ . . 

Z = I — a z a x I cosa;ot — & y sin uj t (21) 

\L0 C LV C J 

which have usual commutation rules: [X, Y]_ = 2iZ,.. In terms of these operators the 
Hamiltonian of the system can be rewritten as 



H = ^X-^Z 
2 2 



A . . e 
— [Z cos uj t + Y sin u t) H X 



Q + H B . (22) 



where A is proportional to the amplitude of the driving force, A = (AF/u> c ). Here we 
resort to the rotating wave approximation (RWA) and neglect fast oscillating terms. Taking 
into account an explicit time dependence of the operators Y and Z we derive the following 
Heisenberg equations (h = 1, k B = 1) 

X = AY + 2— (Y cos Lu t - Z sin iu t) (Q + f), 

Y = 5Z- AX -2 (— Xcoscuot- — z) (Q + f), 

\U C 0J C J 

Z = -5Y + 2(—Xsmuj t- —y) (Q + f). (23) 

In the case of a Gaussian statistics of a free heat bath variables or for a weak qubit- 
bath interaction a response of the heat bath on the action of the qubit is determined by the 
expression [22] 

Q(t) = g (0) (t) + / dtMt,ti)cr z (ti), (24) 



where 



A e 

— [Z cos cu t + Y sin cu t] H X. (25) 
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As in the Section II, a retarded Green function of free heat bath is denoted by tp(t,ti), 
ip(t,ti) = (i[Q(°\t), Q(°\ti)]_)9(t — ti) with a respective Fourier transform (a susceptibility) 
x{w). This susceptibility, xi 00 ) = Xo(^) + Xt(u), incorporates a part Xo( w )) that is due 
to internal mechanisms of qubit decoherence, together with a resonant contribution of the 
tank, xt(uj), 



X t{u) = k?L q I* q 2 7 . ■ (26) 

Besides that, the free heat bath is characterized by a correlation function M(t, ti), M(t, ti) = 
((l/2)[Q(°\t), Q^(ti)] + ), and by a spectrum of equilibrium fluctuations S(cu) with temper- 
ature T, 

S(u>) = J dre^M(r) = X » coth (^) . (27) 

This spectrum S(u) = Sq(u) + St(ou), contains a part originated from the qubit interaction 
with its own heat bath, So(lj), as well as a part, St(uj), related to qubit coupling to the 
tank. According to the fluctuation-dissipation theorem [21], the equilibrium spectrum S(u>) 
is proportional to the imaginary part of the heat bath susceptibility, x"{ u )- 

Following to the method outlined in the Appendix A we can rewrite the Heisenberg 
equations (23) in the form of quantum Langevin equations, 

X = AY + L x + £ c + f x , 
Y = 5Z-AX + L y + i y + f y) 

Z = -5Y + L z + £ z + f z , (28) 

with the collision terms 

2A r 

LJt) = — / dti{M(t,ti)i\Y(t) coso;ot - Z(t) smcu t, 

AZ(ti) coscj ti + AY(ti) sinc^i + eX(ti)]- + 
ip(t, ti)(l/2)[Y(t) cosu t — Z(t) sinu t, 

AZ(*i) cosc^i + AY(ti) sinu^i + eX(ti)]+}, 
2 r 

L y (t) = / dti{M(t,ti)i[AX(t) cos uj Q t - eZ^), 

AZ(ti) coscj ti + AY(ti) sincj ti + eX(ti)]- + 

V?(Mi)(l/2)[AX(f) coscj t - eZ(t), 

AZ(*i) coscj ti + AY(ti) sinc^i + e-X"(ti)]+}, 
2 /• 

L a (i) = — / dti{M(i,*i)i[AX(*)sinw t-eY(t), 

AZ(ti) coscj ti + AF(ti) sinc^i + eX(*i)]_ + 
V?(t,ti)((l/2)[AX(t) sincj t - eY(t), 
AZ(*i) cosw fi + Ay(ti) sincj ti + eX(ti)]+}, (29) 

M(r) = M(t)6(t),t — t — t\, and the fluctuation sources £ x ,£y,£z- Definitions and corre- 
lation functions of these forces, £ m (t) = {Q^°\t), A m (t)}{m = 1,2,3) are presented in the 
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Appendix A (see Eq.(A8)). Hereafter the digital indices 1,2,3 correspond to the indices 
x,y,z, respectively. The qubit operators A m (t) are defined as follows (see also Eqs.(23)): 

A 

A x (t) =2 — (Y cos u t — Z sin ui t), 
u c 

AJt) = -2 (— Xcosuj t - —z] , 

\UJ C UJ C J 

AJt) = 2 (—X sin iu t - —y) . (30) 
\uj c u c J 

We also introduce the effective forces f x ,f y ,fz,'- fm{t) = A m (t)f(t), which are necessary 
for calculating the linear response functions and susceptibilities of the qubit. After the 
calculations the auxiliary force f(t) should be set equal to zero. 

The non-Markovian stochastic equations (28) can be simplified in the approximation of 
weak coupling between the qubit and the heat bath (Bloch-Redfield approximation). With 
non-zero detuning, 5^0, the free evolution of the driven qubit (without coupling to a heat 
bath) is described by the equations (r = t — ti): 

,w s ,w x <5 2 + A 2 cos Q r t . A . ^ r _ r . , A _ 1 — cos il r t 
X(t) = X(h) =j 5_ + Y(ti) — — sm Qrt + ZfaAS 

Y(t) = Y(t x ) cos il R r - XitJ-^-smflRT + Zfa)-^- shifter, 

Mr IIr 

7(+\ vu ^ A 2 + 5 2 cosQrt vU ^S . . 1- cosQrt 

Z(t) = Z[ti) pjj Y(ti)—smQ R r + X[ti)M ^ , (31) 

\L R IIr \l r 

where Qr is an effective Rabi frequency of the qubit, 



(A W \ 2 
—) +(uj -lu c ) 2 . (32) 

In the Bloch-Redfield approximation we can reduce the qubit operators taken at the mo- 
ment t to the operators at the moment t 1 using Eqs. (31), and, thereafter, calculate 
(anti)commutators of the simultaneous qubit operators using usual commutation rules: 
[X(ti),Y(ti)]- = 2iZ(t 1 ),[X(t 1 ),Y(t 1 )] + = 0,.. (see Appendix B). Neglecting fast oscil- 
lating terms in the collision integrals (29) we derive the following equations for the qubit 
operators X\ = X , X 2 = Y, X 3 = Z in the rotating frame of reference (m, n — 1, 2, 3): 

X m + ^2 U mn X n + X! / dtiT mn (t — ti)X n (ti) = + f m + v m , (33) 

n n 

where for nonzero elements of the matrix II we have: II12 = — n 2 i = —A, n 23 = — n 32 = 
—5. The collision coefficients T mn (r) are presented in the Appendix B. For steady-state 
parameters v\ — v x , v 2 — v y , f 3 = v z we obtain the following expressions: 



v x = 2 



A 2 
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1 

2 



5 



1 



^2 



R, 



2d 2 



V 3 



A 



R 



8 A s 

Q R d 2 



x"M+ 



2d 2 



(34) 



Here x'i 00 ) an d x"^) are re& l an d imaginary parts of the heat bath susceptibility xi 00 )- 
A formal solution of the equation (33) has the form 



X m (t) =J2G m n(t)X n (0) + dhG^t - fi)[£„(tl) + fnih) + "n], 



(35) 



where the last constant term describes the steady-state values of the average qubit variables: 
X lfi = X = (5/Q R )P ,X 2 ,o = Y = 0, and X 3i0 = Z = (A/Q R )P , with a polarization 

4(eA/Af x "(Q R ) + (Q R - 5) 2 X "(uj + Q R ) - (Q R + 5) 2 X "(u - Q R ) 

A(eA/Ays(n R ) + (n R - 5) 2 s(d + n R ) + (n R + sys{u - n R ) ' 1 } 

For an exact resonance, c^o — co c ,S — 0, between the frequency of driving force d and 
energy splitting of the qubit d c = \/ A 2 + e 2 , the steady-state polarization P is positive for 
the Ohmic or super-Ohmic heat bath, x"( u ) ~ of ,r > 1. However, P can be negative at 
non-zero detuning, 5^0. 

Here we consider the case of strong driving when the Rabi frequency Q R (32) is much 
more than qubit's relaxation rates. Then, for Fourier transforms G mn (d) of the Green 
functions G mn {r) incorporated into Eq.(35) we obtain: 

GnH = (5 2 -d 2 )/D(d),G 22 (d) = -u 2 /D(u), G 33 (u) = (A 2 - d 2 )/D(d), 
G 12 {uo) = -G 2l {uo) = -iuA/D{u),G l3 {u) = G 31 (cd) = 5A/D(d), 

G 23 (d) = -G 32 (d) = -iu5/D(u), 



(37) 



with a denominator 



D(uj) = i[u + iT z (d)][d 2 -tt 2 R + idT(d)] ~ 
% [uj + iT z (w)] [u - fl R + iT (u) 1 2] [to + VL R + iF(u) /2] . 



(38) 



The coefficients r z (u),r(u), derived with Eqs.(B2),(B3) from the Appendix B, play roles 
of frequency-dependent relaxation rates. These relaxation rates are even functions of 
d,T z (— d) = T z (u), r(— id) = r(u), and they are determined by the spectral density of 
the heat bath S(u) (27): 
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2 A 2 

r,M = — 2 t^\s{uj + 0^) + 5( w - n R )\ + 

2 



(39) 



rM = 2^^) + 



2 X2 r 



2u 2 



(i - -^j s(u + n R ) + (i + -£) 5(cu - 

1-^-) 5(w + w )+fl + ^ 5(w-w ) 



+ 



+ 



^1 + — ") [5 {uj - ujq - Q R ) + S(lu + lu - Q R )] . 
to J 



2ui n R 



(40) 



We omit here frequency shifts of the qubit resulting from its interaction with the heat bath. 
With Eqs.(37) we can calculate the retarded Green functions G mn (r) defined at r > 0: 



G xx (t) = ^e- r - + ^"^cos^r, 



G xy (r) = —G„ x (t 



A 



yx \ 



e" rr/2 snifter, 



R 



G xz (t) = G zx (t) 



5A 



-r z r _ -Ft/2 



COS 



G yy( T ) = e cosQrT, 
G yz (r) = -G zy {r) = A e -rr/2 sinnRT) 



V R 

A 2 A 2 
71 -r z r , „-Tt/2 



G zz (t) = + — e - LT ^cosn R r, 



(41) 



where G mn (0) = 5 mn . Decay rates T z and T are equal to the functions T z {uj) and T{uj) 
(39), (40) taken at zero frequency and at the Rabi frequency, respectively: 



2 A 2 

r, = r,(o) = 2~s(Q R ) + 

u c lL R 



2^ 



R , 



(42) 



r = v{n R ) 



e 2 5 2 



A 2 A 2 
2^^S(cu ). 



(43) 
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The decay rate T/2 is related to the rate T^ 1 from Ref. [18], a notation for Y z remains the 
same. It should be noted also that the heat bath operator, Q(t), defined in the present 
paper differs from the same operator Q from the paper [18], by the factor 1/2 (see also [33]). 
Because of this, to compare our results with results of the above-mentioned article we have 
to divide our spectral function of the heat bath S(u) and our susceptibility x( u ) by four to 
get the spectral function and the susceptibility of the heat bath used in Ref [18]. 



IV. DISSIPATIVE DYNAMICS AND A LINEAR RESPONSE OF THE QUBIT. 



It follows from Eq. (35) that an evolution of the averaged qubit operator (X m (t)),m = 
1,2,3, from its initial condition (X m (0)) to the steady-state value X mt0 is governed by the 
corresponding Green functions G mn (t) (41): 



(X m (t)) = X m , + J2 G mn (t)[(X n (0)} - X nfi ] 



(44) 



In particular, if the qubit starts from the ground state |0) of the Hamiltonian H = (A/2)a x + 
(e/2)a z = (uj c /2)X, where (0|X(0)|0) = -1, (0|F(0)|0) = (0|Z(0)|0) = 0, a probability to 
find the qubit in the excited state, P E — (1 + (X))/2, 



P E (t) = \(l + £p ) (l - |re^ - |^e- r */ 2 cosfi^ 



R 



5A 2 
2Q| 



-r z t 



e- r '/ 2 cos 



(45) 



oscillates with the Rabi frequency Hr (32) and relaxes to the steady-state value PE,st = 
(1/2) [1 + (S/Qr)Po\. In the case of resonant driving when u; — ou c ,5 — 0, the energy levels 
of the qubit are populated equally in the steady state, PE,st — 1/2. At the same initial 
conditions the averaged current in the qubit loop, (I q (t)) = I q (a z (t)), oscillates not only 
with the Rabi frequency Qr, but also with the frequency of driving force c^o, as well as with 
frequencies Uq ± Qr : 



A A 

= — 7T- IqCOSUJ t 



R 



5 



1 + 



-r.t 



e- rt / 2 cos 



n R t) 



+ 



— ——IqSmuJotsmVlRte r<//2 



UJr. 



n 



-p 







R 



R R 



rt / 2 cosft R ^ - 



£ 1 5 p 

LO c ILr 



^( e - r ^- e ^/ 2 cos^t 

R 



(46) 



The qubit starts with an initial current, (I q (0)) = —{e/uj c )I q , corresponding to the ground 
state with a non-zero bias e, and ends at t ^> T~ l ,2/T with the steady-state current oscil- 
lating with the frequency of the driving force Uq : 
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(I q (t)) 



St 



r „ / A A Z $ 

l q P I cosu Q t + 



IjJr 



R 



IjJr ^ 



R , 



(47) 



Interestingly, there are no signs of the Rabi frequency in oscillations of the steady-state 
qubit current. We emphasize that a relaxation of the population difference Pe and a decay 
of the qubit current (I q {t)) are determined by both damping rates T z and r/2 (42), (43). It 
is worth to note also, that zero-frequency fluctuations of the heat bath, such as 1/f- noise, 
which are described by the spectral function S(0), contribute to the decay rate T (43) in the 
case of non-zero detuning S and non-zero bias e. 

Equations (35) averaged over the initial state of the qubit and over the 
thermodynamically-equilibrium initial state of the heat bath allow us to calculate a lin- 
ear response of the qubit on the action of small external force f(t), namely, a functional 
derivative (5a z (t)/5f(t')) or its Fourier transform, a magnetic susceptibility x zz {uj) (10). 
With Eq.(25) we obtain the following expression: 



,5<jJt). A.5Z(t), A.SYCt), . e ,5X(t), 



(48) 



The derivatives of the qubit operators in the rotating frame of reference, {5X m (t)/Sf(t')), 
can be found from the averaged equation (35) taking into account formulas for the forces 
fn(t) '■ f n (t) = A n (t)f(t), where A n (t) are defined by Eqs. (30). The mean values of the 
averaged qubit variables (X), (Y), (Z) in Eqs. (30) should be replaced in the process by their 
steady-state values Xo = (S/Qr)Pq,Yq = 0, and Zq = (A/Q r )Pq, where the polarization 
Pq is given by Eq. (36). Then, for the magnetic susceptibility of the qubit we obtain the 
following result: 



A*_ 

4wr 



1 + 2 



5(uj + ujq 



UJ 



e^A 2 
~u^Q 



1 



R 



4w* 



Po 



1-2 



O 2 
5(lo- 



+ 



i 



1 



Ok 



uj + uj — + iT/2 uj + lo + Q R + iT/2^ 
( 1 1 



uj — ujq — £l R + iT/2 uj — uj + Q R + iT/2 



(49) 



The imaginary part of this susceptibility which defines absorption properties of the driven 
qubit peaks at the Rabi frequency Q R as well as at frequencies uj ± Q R : 



2 A2 
II ( \ ZD 



R 



r/2 



r/2 




4a, 2 



( w _ n R )2 + (r / 2 )2 (c + ^) 2 + (r/2) 2 
r/2 r/2 



uj + lu - n R y + (r/2) 2 (lu-lu + n R y + (r/2) 2 
r/2 r/2 



+ 



+ 



R, 



(uj-uj q - n R f + (r/2) 2 {uj + uj + n R y + (r/2)< 



(50) 



An absorption of weak signal energy by the qubit is determined by the function U{ui) = 
ujx zz (uj) [34]. This is evident from the formula (50) that U(u>) can be negative at the positive 
frequency uj — Q R . It means that a weak signal having this frequency will be amplified by 
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the strongly driven qubit. The low- frequency part of the qubit susceptibility which affects 
the resonant frequency of the tank (see Eq.(19)) has the form 



xLM = - 



UJ~Q R 



UJ + fl R 



(lu - n R y + (r/2)2 ( w + n R y + (t/2)< 



(51) 



The angle 9 between a voltage and a current in the tank coupled to the qubit (14) is 
determined by the susceptibility of the qubit (50), (51) taken at the tank frequency ojt '■ 



tan0 = -k 2 L n I. 



,e 2 A 2 



-Pi 



R — UJT 



q "u 2 n 2 R \n R -u T y + (T/2f 



It' 



Here jt, 



e 2 A 2 

7t = It + k 2 L q I 2 — -^P uj T] 



r/2 



(n R - uo T f + (r/2)2 

is an effective damping rate of the tank in the presence of the qubit. 



(52) 



(53) 



V. NONEQUILIBRIUM SPECTRA OF THE QUBIT AND THE TANK. 

Here we will calculate the total nonequilibrium spectrum of the qubit fluctuations, S zz (u>) 
(17), for the case of exact resonance between energy splitting of the qubit, uj c = \/ A 2 + e 2 , 
and the frequency u of the driving field: 5 = u — u c = 0. Calculations of the spectrum 
for non-zero detuning are straightforward, but cumbersome enough. The part of the qubit 
spectrum S ZZ:0 (u), which results from internal decoherence mechanisms of the qubit (cou- 
pling to the bath Q ), can be easily found from the expression for the total spectrum S zz . 
To do that we have to replace the total spectrum of the bath, S(lu) by the spectrum of the 
internal bath Sq(u) in the expressions (C2)-(C7) for the spectra of fluctuation forces given 
in the Appendix C. It should be emphasized that damping rates of the qubit are determined 
nevertheless by the total spectrum S(u) (27) of the dissipative environment. 

In view of Eq.(25) the correlator of cr 2 -operators of the qubit averaged over fast os- 
cillations can be expressed in terms of qubit's correlation functions in the rotating frame 
(r = t-t>): 

(l[a z (t),a z (t%) = ^[X(t),X(t%) + 

A " ^ [Z(t), Z{t')] + ) COSLUoT + (1 [Y(t), Y(t>)} + ) COS U!qT-\- 

{ l - [Y(t), Z(t')} + ) sin^r - (I [Z(t), Y(t')} + ) sin^ r} . (54) 



Correlators of the qubit variables in the rotating frame are determined by the correlation 
functions of the fluctuation forces (see a stochastic part of Eq.(35)) 

(i [X m (t),X k (t')} + ) = J ^AE^(Ml)G H (M2)(J [£n(*lUi(*2)]+), (55) 
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where the Green functions, C mn (r), are defined by Eqs.(41). Fourier transforms of the qubit 
correlation functions, A mk (uj), and the correlator of the fluctuation forces, K n i(u), 



(\[X m (t),X k (t)] + ) = J ^ e -M*-0 Amfe(w); 



{- 2 [Ut)Mt')\ + ) = J^e~^K nl {ul 
are related according to the equation 

A rnk (u) = ^2 G mn (u)G k i{-u)K n i(u) (56) 

lq 

with G rnn (uo) from Eqs.(37). Correlators of fluctuation forces, ((1/2) [£n(t),£i{t')] + ), are cal- 
culated according to the procedure given in the Appendix A (see Eq.(All)). The expressions 
for spectral functions of the fluctuation forces, K n i(u>), are presented in the Appendix C. 
For the spectrum of qubit fluctuations S zz (uj) we find from Eqs.(17),(54): 

e 2 

S zz (uj) = —A xx (u) + 

A 2 



4w* 
A 2 



4cu 2 



[A zz (u - u ) + A yy (uj - uj ) + iA yz (u - u ) - iA zy (u + u )]. (57) 



It follows from Eqs.(57),(37), that the spectral functions of the qubit operators in the rotating 
frame, A mk {uj), (m, k — x,y, z) are determined by the spectra of fluctuation forces K n i{uj): 

A xx (lu) = lu 2 [lu 2 K xx (lu) + n 2 R K yy (iu) - 2iuAK xy (tu)]/\D(Lu)\ 2 , 
A yy (uj) = u 2 [n 2 R K xx (u) + u 2 K m (u) - 2iu AK xy (u ) ] / 1 D (u ) | 2 , 

A zz (co) = (co 2 -n 2 R ) 2 K zz (uj)/\D(uj)\ 2 , 
A yz (u) = iu(Q 2 R -u 2 )[AK xz {u)+iu;K yz (u;)}/\D(u;)\ 2 , (58) 

with A zy (uj) = —Ay Z {uj). Here 

|/JH| 2 = [uj 2 + r 2 M][(c 2 - n 2 R ) 2 + c 2 r 2 M] (59) 

is the modulus square of the Green function denominator (38). Combining Eqs.(57)-(59) 
with the formulas (C2)-(C7) from the Appendix C we obtain the nonequilibrium spectrum 
of qubit fluctuations S zz (cu) : 

s m g2 Wr{uj) I 

zz[ ' 2lu 2 (lu 2 - n 2 R ) 2 + lu 2 T 2 (lu) + 
A 2 W(uj + u) 1 

X 77 77 „o m 77-77T7 7 



4u 2 (u + u ) 2 + Tl{u + uj ) [{u + u ) 2 - Q 2 R } 2 + u> 2 r 2 (u + u ) 

A2 W ^ - g) x I (60) 

4uj 2 (u-u ) 2 + r 2 (u-u; ) [{io-uj ) 2 -tt 2 R } 2 + u 2 T 2 (u-Uv) 1 ' 
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Here the frequency-dependent damping rates T z (uj),r(uj) are defined by Eqs. (39), (40). 
Functions W R {ui) and W{uj) are given by the expressions 

W R (u) = 8Q|4^H + ^ 2 ^{S{oo + cu ) + S{uo - coo)} + 

^ c ^ c 

A 2 

— (u - Q r ) 2 [S(uj + Uq + Q r ) + S(u - uj + Qr) - 



P X"(U + LOQ + ton) - P0X(U -UJQ + Ofl)] + 

A 2 

— (w + n R ) 2 [S(cu -co - Q R ) + S(cu + cu - n R ) + 



P oX "(lu -lu - Q r ) + P oX "(lu + luo- n R )], (61) 



W(w) = Auo^S{uo) + uo 2 n R ^[S{uo + cu ) + S{uo - UO Q )} + 



£ 2 



i—uj - q 2 r ) 2 [S(uj + - PoxV + + s(uj - n R ) + p x> - n R )] + 



c c 2 



A 2 

^r^(^ - fi R ) 2 [.s(^ + cu + - p oX "(uj + u + n R )} + 

A 2 

Q2r 2^ ^ + + ^° ~ ^ + + ^ _ ^ + 

A 2 

— (a; - R ) 2 (2a; + fi^) 2 ^ - cu + - Pqx'V - + ^)] + 

A (^ + fi i? ) 2 (2^-fi i? ) 2 [ 1 S(^-cu -fi R )+PoX ,/ (^-^o-fi R )]. (62) 



2a; 2 

It follows herefrom that the qubit spectrum is an even function of uj : S zz (—ou) = S zz {uj). 
In the regime of strong driving, when Q R ^> r,T 2 , the qubit spectrum S zz (u) taken at 
positive frequencies consists of three Lorentzian peaks centered at frequencies Q R , oj q — fl R , 
and Uq + VL R with the linewidth T/2 each, with T = T(Q R ) (43), complemented by the 
additional peak located exactly at the frequency of the driving force c^o- In the present 
section ojq is equal to the frequency of quantum beatings u c . The additional peak has a 
different linewidth, T z = T 2 (0) (42), and it is absent in the absorption spectrum determined 
by the function x'Li 00 ) (50). The current in the qubit loop is described by the operator 
I g = I q cr z . Because of this the spectrum of current fluctuations of the qubit, Si(cu), is 
proportional to the spectrum S zz {uj) : Sr(ou) = ((l/2)[I q (u), I q ]+) = I 2 S zz (u), where I q is 
the value of the persistent current in the qubit loop. 

In view of the facts, that cu > Q R , and W R (Q R ) ~ 4TQ 2 R , W(0) ~ 2T z tt R , W(Q R ) ~ 
2TVL R , and considering the heat bath with temperature that is greater than the energy of a 
Rabi quantum, T ^> hQ R , we find the simple formula for the nonequilibrium spectrum of 
qubit fluctuations: 

^_ n e 2 n 2 R T 
bzz[U} ~ uU^ 2 -^ 2 r ) 2 + uj^ + 
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A 2 T z A 2 r 2 

%4 {uo + coo) 2 + r 2 + 2^f(cu-^ ) 2 + r 2 + 
r r 

2ul [{uj + cu ) 2 - n 2 R ] 2 + cu 2 P 2cu 2 [(u - uoo) 2 - ^l? + ^ 2 P ' 1 J 

where the decay rates can be found from Eqs.(42),(43) at zero detuning 5 = 0: 

r z = 2^s(n R ) + ^s(u ), 

r = 2^S(n R ) + 3%S(u ). (64) 

As for the driven atom [16] the spectrum of our two-level system is double peaked at the 
frequencies o;o ± VL R , but, besides that, we have a peak at the Rabi frequency with the 
intensity that is proportional to the bias squared e 2 . The low- frequency part of the spectrum 
S zz (uo) gives a significant contribution to the voltage spectrum of the tank Sy(u) (18). We 
recall also that the internal heat bath only should be taken into account in the process of 
calculating the spectra of fluctuation forces K n i(uj) (C2)-(C7). It means that in expressions 
(61), (62) for W R (u), W(uj) we have to extract fluctuations of the tank from the total heat 
bath spectrum, S(u>), and substitute this spectrum for the spectrum So(u) related to the 
internal bath. The damping rates V and T z in numerators of Eq. (63) are originated exactly 
from the function W(u) (62). For calculating the spectrum S ZZj o(iv) it is necessary to replace 
these rates by the coefficients depending on the spectrum of the internal bath only, namely, 
r 2i0 and r , where, for example, T = 2(e 2 /uo 2 c )Sq(VL r ) + 3(A 2 /o; 2 )S , o(cc;o). The decay rates 
in the denominators of the spectrum S ZZ fi{ui) remain the same, because both the internal 
mechanisms and the tank fluctuations contribute to the linewidth of the qubit. As a result, 
for the low-frequency part of the qubit spectrum S zz p{uj) we find 



It should be noted that all of these nuances with replacing S(u) to So (a>) and T to T in 
the spectrum S zz (cu) (63) are important only when the contribution of the tank into the 
fluctuations and decoherence of the qubit is quite significant. For weak qub it-tank inductive 
coupling a strong influence of the tank on the qubit coherence can take place only near the 
exact resonance between the Rabi frequency Q R and the resonance frequency of the tank 
ujt- It is necessary to develop more rigorous theory to study a close proximity of this point. 

The total spectrum of voltage fluctuations in the tank, Sy(uj) (18), incorporates a contri- 
bution of the tank noise, Svt(w), together with a low-frequency contribution of the driven 
qubit, Svq{w) : Sy(u) = Syri^) + Svq{u). At temperatures T ^> uj^ the tank contribution 
to the voltage spectrum is proportional to the spectral function S(,(iv) = T^ t /ujt (7): 



, u 2 T 7T 
fe(^-cu 2 ) 2 + ^ 2 7 2: 



Syriu) = 2—— ^ —. (66) 



where ojt is the resonance frequency of the tank (19) shifted due to qubit-tank coupling, 
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UT = UT f + PL ^5 P ^c, T -X) 2 +V/2)^ (67) 

We use here Eq.(51) for the function x'(u;t) -The linewidth of the tank, 7^, modified by the 
qubit, is defined therewith by Eq.(53). 

For the qubit part we find from Eqs. (18), (65): 

e 2 L I 2 uj 2 f2 2 

Svq(v) = -p^-UJ r 0T-T2 J X T^2 n2 \2 i 2^2' ( 68 ) 

cj 2 (ui t - cj 2 ) 2 + w 2 7^ (cu 2 - f2jj 2 + a; 2 r 2 

Measurements of the voltage fluctuations are performed within the linewidth of the tank: 
uj ~ Co T ± 7 T . In this frequency range a signal-to-noise ratio demonstrates a resonant be- 
haviour as a function of the Rabi frequency Vl R : 

S VT (u) ]ui=U)t T lT (u 2 -Q 2 R ) 2 + uj 2 r 2 - { > 

In this expression we have a ratio of bare damping rates of the qubit, T , and the tank,7^. 
Besides the part To, related to the contribution of the internal heat bath into qubit deco- 
herence, the total decay rate of the qubit T (64) contains also a term, IV, which describes 
a tank share in the qubit damping: r = T + T T , with 

T T = Ak 2 L q I 2 q — 2 uo 2 T 2 -^g— - 2 2 . (70) 

This rate reflects a backaction of the detector (LC circuit) on the quantum bit that accom- 
panies an acquisition of any information from the qubit. Both parameters, Svq/Svt and 
IV reach the maximums when the Rabi frequency Qr is about the resonant frequency of 
the tank u T . However, the signal-to-noise ratio (69) as a function of the Rabi frequency 
VLr has a linewidth that is determined by the damping rate of the qubit T, whereas the 
tank contribution to the qubit decoherence (70) is localized in the narrower range of Rabi 
frequencies which is of order of the tank linewidth 7^, 7t "C V. Measurements performed 
at the Rabi frequencies that are out of this range, VLr > ojt + 7t, or VLr < ujt — 7t, 
can demonstrate a good efficiency, Svq/Svt > 1, without introducing strong decoher- 
ence in the qubit. To observe a signature of Rabi oscillations in the spectrum of volt- 
age fluctuations Sy(ou) the decay rate T and, all the more, the measurement-induced rate 
IV should be appreciably less than the Rabi frequency Qr ~ ujt '■ T^/out <C 1- For the 
flux qubit measured in Ref. [7] we have the following set of parameters: L q = 24pH, I q ~ 
500nA,Q T = uj t /2 1t = 1850, T = lQmK, lu t /2it = 6.28AM Hz,hu T = 4.16 x 10- 27 J, so 
that Tjhw T = 33, L q I 2 = 6 x 10~ 24 J, L q I 2 /hu T = 1440, co T /-f T = 3700. If we take a value of 
the coupling parameter squared, k 2 ~ 10~ 3 , from Ref. [7] and suppose that e/A~ 1/126, 
then for the measurement-induced damping rate we obtain the ratio: Tt/ujt — 0.8 x 10 -2 , 
at the point where VLr — ujt — T/2 ~ O.Olu;^. We use here the decay rate, T = 0.02u;t <C wt, 
measured in Ref. [7]. At these conditions the signal-to-noise ratio (69) is of order 0.5, and 
the detector-induced decoherence of the qubit, T T , as well as the total rate T are much less 
than the Rabi frequency of the qubit Qr. It means that the spectroscopic observation of 
Rabi oscillations [7] can be classified as a weak continuous quantum measurement. 
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VI. CONCLUSIONS. 

In this paper we have analyzed quantitatively a continuous spectroscopic measurement 
of Rabi oscillations in a flux qubit by means of a high quality tank (LC circuit) which is in- 
ductively coupled to the qubit loop. This circuit serves as a linear detector for measuring the 
spectrum of voltage fluctuations in the tank as well as for monitoring the effective impedance 
of the system " qubit +tank" . Expressions for the voltage spectrum and for the angle be- 
tween the current that drives the tank and the averaged tank voltage have been derived in 
terms of the spectrum of qubit fluctuations and the qubit magnetic susceptibility. To find 
the spectrum of the qubit and its magnetic response we have applied a formalism of non- 
Markovian Heisenberg-Langevin equations to the case of a strongly driven open quantum 
system. Combining the Bloch-Redfield and rotating wave approximations we have obtained 
formulas for the damping rates of the qubit and its magnetic susceptibility as functions of 
the amplitude of the high-frequency driving field and detuning of this field from the qubit 
energy splitting. A dissipative evolution of the averaged current in the qubit loop and the 
probability to find the driven qubit in the excited state have been described analytically. 
Contributions of the qubit to the damping rate and the frequency shift of the tank have been 
calculated. We have presented also analytical formulas for the nonequilibrium spectrum of 
current fluctuations in the qubit loop as well as for the spectrum of voltage fluctuations in 
the tank (detector) which contains information about the Rabi frequency and about the 
decay rate of Rabi oscillations T. It is shown that the ratio between the qubit contribution 
to the spectrum of voltage fluctuations and the thermal spectrum of the tank is peaked 
when the Rabi frequency is about the resonant frequency of the tank lot- It corresponds 
to the maximal acquisition of information from the qubit. We have shown also that this 
effective measurement is accompanied by the maximal value of decoherence resulting from 
the backaction of the tank on the qubit. The signal-to-noise ratio as a function of a deviation 
between the Rabi frequency and the frequency of the tank has a linewidth that is propor- 
tional to the qubit decay rate T, whereas measurement-induced decoherence of the qubit 
as a function of the same deviation is determined by the linewidth of the tank 7^ which 
is much less than T. It allows us to find optimal conditions for the efficient spectroscopic 
measurement of the Rabi oscillations in the strongly driven quantum bit. 
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APPENDIX A: METHOD OF QUANTUM LANGEVIN EQUATIONS. 

In this Appendix we sketch out basics of our approach to the theory of open quantum 
systems which has been proposed in Ref. [22] and developed in Refs. [23,24]. Formulas for 
the collision terms L x ,L y , L z (29) will be derived in the process. Besides that, we give here 
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explicit expressions for the fluctuation sources £ x ,£y, £z together with a recipe for calculating 
their correlation functions. 

The Heisenberg equations (23) incorporate the total heat bath operator Q(t) (24) multi- 
plied by an operator of the qubit, say, A(t). These operators commute because they belong 
to the different physical systems. It is convenient to work with the symmetrized product of 
these operators. With the expansion (24) in mind we obtain 

\[Q(t),A(t)] + = \[Q< 0) (t),A(t)] + + J dtMtM)\[° z {h),A{t)\ + . (Al) 

The averaged value of the first parametric term in this expression is determined by the 
quantum Furutsu-Novikov theorem [22]: 

(\[Q^{t\A{t)] + ) = J dt'M(t,t')(^L), (A2) 

where M(t, t') is the symmetrized correlator of unperturbed heat bath variables, M(t, t') = 
((l/2)[Q(°)(i),Q(°) (?)]+). This relation follows from the fact that due to the qubit-bath 
interaction a Heisenberg operator of the qubit represents a functional of the bath variables 
{Q^} ■ A(t) = A[{Q^},(t)] that can be expanded in a sum of various time-ordered 
products like Q {0) (t ai )Q^(t a2 )...Q^ (t a J. The operators of the free heat bath, Q(°), obey 
the Wick theorem. Because of this the average value of the additional operator Q^(t) 
multiplied by the term Q^(t ai )Q^ (t a , 2 )...Q^ '(t an ) is reduced to the sum of pairings between 
the external operator Q^(t) and each operators from the above-mentioned product: 

(Q m (t).Q^(t ai )Q^(t a2 )...Q^(t an )} = (Q {0) (t)Q^(t ai )}(Q^(t a2 )...Q^(t a J}+ 

(g(°)(t)g( )(t a2 ))(g( )(t ai )g( )(t a3 )...g( )(t a j) + ... 

The operator (t ak ), k = 1, .., n, engaged with the operator (t) in the outside correlator 
(Q^ (t)Q(°) (t ak )) , should be removed from the initial product. It corresponds to taking a 
functional derivative over this variable: 

(Q(°\t).Q( \t ai )Q( \t a2 )...Q(°\t an )) = 

/^(g (0 H0Q (0 HO>(^|^{Q (0) (^)g (0) (ia 2 )...g (0) (iaj}>. 

This equation together with a relation 5Q^(t a )/5Q^(t') = S(t a —t') results in the formula: 
(Q^(t)A(t)) = J ^(Q(o)( t )Q(o) (0)( _^|_ > . (A3) 

We notice that the position of the external operator Q^(t) with respect to the operator A(t) 
is mapped into a relative order of operators in the commutator (Q^ (t)Q^ (£')) involved in 
Eq. (A3); in so doing the average value of the symmetrized product of Q^(t) and A(t) (see 
Eq.(Al)) is determined by the symmetrized correlator of the heat bath M(t,t'). 
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The functional derivative over the variable Q^(t') is equivalent to the derivative over 
the deterministic force fit') which is additive to Q^(t') in the Hamiltonian (1). In its turn, 
the functional derivative of the qubit operator A(t) over the force fit') is proportional to 
the commutator of A(t) and the qubit matrix cr z (t') that is conjugated to the force fit') in 
Eq.(l): 

SA(t) _ SA(t) _ . [A{tl az{tl)] _ e{t _ t , y (A4) 



6QW(t) 5 f{t') 

To show this we consider the Heisenberg operator of the qubit, Ait), in the interaction 
representation, when the interaction between the qubit and the force fit) is described by 
the term H int = —a z f(t). Then, an evolution of the operator Ait), 

Ait) = S + it)A^it)Sit), (A5) 

is determined by the S-matrix, Sit) = Sit, — oo), with 

S(t,t ) = T{exp[* f Wf(fi)/(ti)]}. (A6) 



Here T is a time-ordering operator, af\t), A^\t) are the free qubit operators (without the 
interaction with the force fit) ). Then, for the functional derivative we obtain 

&4(t) SS + (t) {0) s+(t)A Q) (t) MW 
5 fit') + 6 [t) 5fit'Y 

where 

^ = i9(t - t')T{af\t') exp[. dt^O/ft)]} = 
i0(t-t')T{a?\t')exp[i f d*i^ 0) (*i)/(ti)] exp[* f dha^ (ti)/(fi)]} = 



i9(t - t')T{exrt / dha^itMh^it'meMi / ^i^ 0) (ti)/(ti)]} = 

Jt' J-oo 

i9(t - t')S(t, t')a^(t')S(t') = i9(t - t')S(t)S + (t')a^(t')S(t') = i9(t - t')S(t)a z (t'), 

and o" 2 (t') is the total Heisenberg operator. Here we apply the relation S(t,t') = 
S(t)S-\f) = Sit)S + it'), which follows from the facts that S' 1 ^) = S+(t'), and 

Sit, t')S(t') = Sit, t')S(t', -oo) = Sit, -oo) = Sit). 

Taking into account the derivative of the matrix S + : SS + it)/Sfit') = —i9(t — t')a z (t')S + (t) , 
we obtain Eq.(A4) for the functional derivative of an arbitrary Heisenberg operator Ait). 

In view of the Furutsu-Novikov theorem (A2) the operators like (l/2)[Q(t),A(t)] + in- 
volved in the Heisenberg equations (23) can be splitted into a fluctuation force £4 and a 
collision term : 

(l/2)[Q(t),A(t)] + = U + L A , (A7) 
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where the fluctuation force, 

Ut) = {Q (0 \t),A(t)} = \[Q^(t),A(t)} + - J dt'M(t,t') 5 ^, (A8) 
has a zero average value, (£4) = 0, and the collision term, 

L A (t) = J dt ± M(t, ti)i[A(t), a z (ti)\- + J dt^{tM)\[A{t), <r z (*i)]+, (A9) 

incorporates contributions both parametric fluctuations and a back action of the heat bath. 
Here we introduce a causal correlation function of the free heat bath M(t, t x ) = M(t, ti)9(t — 
ti) having S(u>) as a Fourier transform, 

S(u) = I dre^M(r) = f ^ l -—-S{u 1 ), (AlO) 

with S{uj) being the equilibrium spectrum of the heat bath (27), and e — > +0. 

The explicit form of the fluctuation sources allows us to find their correlation functions. 
To do that we have to take pairings of all free heat bath variables with the heat bath 
variables and the qubit operators belonging to other fluctuation forces. In the case of weak 
qubit-bath coupling we can take into account pairings between the free heat bath variables 
only. With this procedure we derive the following expressions for a correlator of fluctuation 
forces U(t) and f B (i): (&i(t)f B (0) = (Q (0) W> Q {0) (t')){A{t)B(t')), and for the symmetrized 
correlation function: 

{\ [Ut),Ut')} + ) = M(t,t')(\ [A(t),B(t')} + ) + R{t,t'){\ [A(t),B(t')}_) (All) 

where the antisymmetrized correlator of the heat bath, R(t,t') = ((l/2)[Q(°)(t), 
has the spectral function x"( w ) as its Fourier transform. 



APPENDIX B: COLLISION INTEGRALS. 

Collision terms L x ,L y ,L z (29) can be simplified in the approximation of weak coupling 
between the qubit and the heat bath. In this case (anti) commutators of the qubit variables 
X, Y, Z taken at different moments of time are calculated with free evolution operators of the 
qubit (31). Here we present expressions for (anti) commutators of qubit operators X, Y, Z 
involved both into the collision terms and into the correlation functions of the fluctuation 
forces. With Eqs. (31) and the usual commutation rules we obtain ( here t — t — t-y): 

A 8 A 

i[X(t),X(ti)]_ = 2Z(t 1 )— smn R T-2Y(t l )- w (l -cos^r), 

1 5 2 A 2 

-[X(t),X(tx)] + = -j - + — C osQ R r; 

Ft Ft 

( 8 2 A 2 \ 8 A 

i[X(t), y(tx)]_ = -2Z(h) f _ + _ cosier + 2X(t 1 )-p(l - cos Q R r), 

\ Ft Ft / Ft 
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-[X(t),Y( tl )l 



l[Y(t),X(t 1 )] + = ^-smn R r; 



i[Y(t), X(h)}- = 2Z(h) cos Q r t - 2Y(t 1 ) — sm Q R r; 

-2 - + -j- cos fi R r - 2X(t x ) — sin Q R r, 

R R / R 

kxitlZ^U = i[Z(f),X(fi)]+ = ^(1 -cos^r), 



^1 



A 2 



= -2y(*x) I ^r + ^rCosfi R r) -2Z{t 1 ) — smQ R r; 



s 



R R / 

i[Y(t), Y(h)}- = 2Z{t 1 )—smtt R T + 2X(t 1 )-^-smn R r, 

\i R \t R 

^[Y(t),Y(t 1 )] + = cosn R r; 

i[Y(t), Z(ti)\- = -2X(t 1 )cosQ R r-2Y(t 1 )—smn R T, 

Mr 

±[Y(t),Z( tl )] + = ~[Z(t),Y( tl )] + = ^smQ R r; 



i[Z(t),Y(t 1 )]_ = 2X(t 1 ) 



A 2 5 2 



SA 



i[Z(t), Z(tx)]_ = 2X(t!)-^ sin^r + 2Y(t 1 )^(l - cosQ R r), 

R R 

[A 2 5 2 n \ 



-[Z(t),Z( tl )}. 



(Bl) 



Using these formulas we find the following expressions for the collision coefficients in- 
volved in Eq.(33): 

5 2 \ „ 5 . . „ A 2 

o 



T xx (t)=2^M(t) 



1 + —y cos fl R r cos ujqt + 2—— sin wot sin f2#T + — -y cos c<JoT 



f x ,(r) = -f ^(r) = 2 



R / 

A 2 A 



6 i? 



cj 2 fi. 



M(r) 



f^(r) = T 2X (t) =2^^-M{t) 
uj \l R 



sin Q R r cos cj t + 7— ( 1 — cos Q r t) sin UoT 
sin r t sin Wot — 7— ( 1 — cos Q r t ) cos c<j r 



£ 2 / A 



u 2 \ n 2 m 



T yy (r) = 2M(r 
f yz (r) = -F zy (r) = 2M(r) 

T zz (t) = 2M(r 



A 2 



<5 2 A 2 



COS CUqT 



COS £l R T 



+ TZ2 COS tt R T 

5 2 5 . „ A 2 / <5 2 A 2 N 

2——— smiZ^r H smuj r —r + —y cos \l R r 



uj 2 tt R 
2 

2 



R R 



c 

A 2 / <5 2 A 2 N 

2— cos iZ^r H coscj r — y + — y cos \l R r 



(-Of. 



^"^ R R 



(B2) 



Fourier transforms of the collision coefficients T mn (r) (B2), 

f m „(u;) = J dre^T mn (r), 
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are proportional to the causal spectrum of the heat bath S(u) (A10). Frequency-dependent 
relaxation rates r z (u),r(u) (39), (40) defining relaxation and decoherence of the qubit are 
expressed in terms of real and imaginary parts of the functions T mn (u) : 

5 2 5A~ A 2 ~ 

Ft Ft Ft 

A 2 ~ ~ 5 2 ~ 

J) A _ & A ~ 

2-f»» + 2-f»» - 2° w Y' xz {u). (B3) 

Lu Lu 1 L R 



APPENDIX C: CORRELATORS OF FLUCTUATION FORCES. 

In this appendix we adduce formulas for the spectral functions of fluctuation forces K n i(u) 
that eventually determine the nonequilibrium spectrum of qubit fluctuations S zz (u) (57). 
Correlation functions of fluctuation forces, £ m (£) = {Q^ (t) , Am(t)} , are calculated according 
to Eq.(A9) with (anti) commutators presented by Eqs.(Bl). For the spectrum K yy {uj), as 
an example, the corresponding correlator of the fluctuation forces, ((1/2) [£ y {t),£ y (t')] + ) is 
obtained from Eq.(A9) with the operators A = A y (t), B(t') = A y (t') (30): 

2^ [ M (t, t'){\ [X(t),X(t')} + ) + R(t, t')(\ [X(t), X(f )]_>} cos u, (t - if) + 

4 5 { MitX){ \ [Z{th Z{t ' )]+) + RM{ \ Z{t ' )] - ) } ■ (C1) 

Following to this procedure for all correlation functions of the fluctuation forces and taking 
corresponding Fourier transforms we find for the spectral functions of fluctuation forces 
K lq (u);l,q = x,y,z: 

A 2 

K xx {uj) = -z[S(u + u ) + S(u - cu )} + 



u 2 c 
A 2 

— — [S(u + lj + tt R ) + S(u - lj - Q R ) + 
S(uo + lj - Q R ) + S(lj - LU + Qr)] + 

A 2 

—P [ x "(co -u - n R ) - x\u + + Qr) + 



2u 



c 



X "(lo + lo - n R ) - X "(oj -lu + Q R )\; (C2) 



A 2 A 

K xy (io) = -K yx (uj) = -i—r — [S(uj + uj + Q R ) - S(u-u - Qr) + 

ZU C \Lr 

S(lJ - Ld + Q R ) - S(UJ + LJ - VLr)} + 
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A 2 A 
2w 2 il R 

X "(uj + uj - n R ) + X "(uj -u + Q R )]; (C3) 



A 2 A 

5(^ + ^0 + ^r) - 5(0; - Uq — Q R )] + 

A 2 A 

2^Q^ P o\x"(v + wo - fifl) - x'V - ^o + fifl) + 

X> + coo + ft*) - x> - w - fifl)]; (C4) 



s 2 „, , A 2 

S(u + ujq - + 5(w - cj + n R )] + 

A 2 

— P [x"(uj - w - - xV + ^o + fifl) + 

x > + w - - X > - ^o + fifl)]; (C5) 



A 2 

K yz (u) = K zy (u) = -i—[S(uj + uj + Q R ) - S(u -u Q - Vl r ) + 

s(uj + ujq - n R ) - s(uj - lo + n R )\ - 

A 2 

i^P G [x"{oj + u - n R ) + x"{oj - uj + n R ) - 

X "(u + u + Q R ) - X > -u - n R )}; (C6) 



e 2 



K zz (u) = 2—[S{uj + + S(u - Q R ) - 
w 2 

p x"(u + n R ) + p oX "(uj-n R )] + 

A 2 

^^[^(w + ^o + CLr) + S(uj - u - Q r ) + 
+ u - tt R ) + S(u -0J0 + n R )] + 

A 2 

—P [x"(u -u - n R ) - x"{oj + u + tt R ) + 

X "(uj + uj - n R ) - X "(uj -u + Q R )]. (C7) 
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